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We report various many-body theoretical approaches to the nonlinear decay rate and energy 
loss of charged particles moving in an interacting free electron gas. These include perturbative 
formulations of the scattering matrix, the self-energy, and the induced electron density. Explicit 
expressions for these quantities are obtained, with inclusion of exchange and correlation effects. 



I. INTRODUCTION 

The energy loss of non-relativistic charged particles entering a metal is primarily due to the creation of electron- hole 
pairs and collective excitations in the solid, interactions with the nuclei only becoming important when the velocity 
of the projectile is much smaller than the mean speed of the electrons in the solidtl. 

The degenerate interacting free electron gas (FEG) provides a good model to describe a regime in which electrons 
are responsible for the energy-loss process. The inelastic decay rate and energy loss of charged particles in a FEG 
have been calculated for many years in the first-Born approximation or, equivalently, within linear response theory. 
It is well known that these first-order calculations predict an energy-loss that grows with the square of the projectile 
charge, Zie, and provide a good approximation when the velocity of the projectile is much larger that the average 
velocity of the target electrons. However, when the velocity of the projectile decreases non-linearities become apparient. 
An important example is provided by the existing differences between the energy loss of protons and antiprotonsBc, 
which cannot be accounted for within linear-response theory. These differences were then successfully accounted on 
the basis of second-order perturbative calculations that used the random-phase japproximation (RPA) and treated 
the moving charged particle as a prescribed source of energy and momentum. trEl Beiiond-RFA calculations of this 
so-called effect have been reported only very recently, in the limit of low velocitiesJlj 

In this paper, we report various many-body theoretical approaches to the quadratic decay rate and energy loss of 
charged particles moving in an interacting FEG, which include exchange-correlation (xc) effects and treat the moving 
charged particle as part of the many-body interacting system. First of all, we present a fully quantum treatment of 
the probe particle, which we assume to be distinguishable from the electrons in the Fermi gas. We assign a propagator 
to this particle, and then follow procedures of many-body perturbation theory to derive explicit expressions for the 
scattering matrix. From the knowledge of this matrix, both the decay rate and the energy loss of the moving particle 
can be evaluated either within the RPA or by including short-range xc effects. We also derive an explicit expression 
for the self-energy of the probe particle, which enables us to present an alternative derivation of the decay rate. 

The rest of this paper is organized as follows: A diagrammatic analysis of the decay rate of a moving charged 
particle in a FEG is presented in Section II. The decay rate and stopping power are calculated up to third order in 
the projectile charge from the knowledge of the scattering matrix. It is shown that for a heavy projectile the decay 
rate agrees with the imaginary part of the projectile self-energy, and that the the stopping power agrees with the 
result of quadratic-response theory. Our conclusions are presented in Section III. Atomic units are used throughout, 
i.e., = ?i = me = 1. 

II. DIAGRAMMATIC ANALYSIS 

We consider the interaction of a moving probe particle of charge Zi and mass M with a FEG of density n. The 
probe particle is assumed to be distinguishable from the electrons in the Fermi gas, which is described by an isotropic 
homogeneous assembly of interacting electrons immersed in a uniform background of positive charge and volume V . 
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In the representation of second quantization, the interaction-picture perturbing Hamiltonian reads 



H'j{t)^-Zi I (frd'^X'ilj''{x)ip{x)v{x,X)'iljHx)ij{X) 

d^rd*x' ij\x) V^a;) v{x, x') i^^x') ^pix') + Hf^, 



(2.1) 



where v(x, x') [x — (r, t)] is the instantaneous Coulomb interaction and the last term represents the interaction of 
electrons and probe particle with the positive background. The field operators ^p{x) and ip''{x) destroy and create an 
electron at time t and point r, while V'(^) and '4'^{X) destroy and create the probe particle at time t and point R. 
Annihilation operators can be written as 



V'(x) =^e^'^'*</.,(r) a, 



(2.2) 



and 



i uJi t 



,(R)A, 



(2.3) 



where the operators Ui and Ai annihilate an electron and the probe particle in the one-particle free states 0i(r) and 
^i(R) of energy uii. As we are dealing with a homogeneous system, these states can be taken to be plane- wave states. 
We choose states of momentum k and energy = k^/2 for electrons, and momentum p and energy ujp ~ p^/(2M) 
for the probe particle. ■— ■ 
The scattering matrix can be written as a time-ordered exponentialllil 



5 = T <^ exp 



I [ dte-"^^'^ H'j{t) \ 

J — OO J J 



(2.4) 



where Hj is the perturbing Hamiltonian of Eq. (2.1), T is the chronological operator, and yy is a positive infinitesimal. 
In the interaction picture, electron and probe-particle propagators can be expressed as 



and 



G{x,x') 



D{X,X') 



. (0,$o|TV>/(a;)v4(a;05|0,$o) 
(0,$o|^|0,$o) 

(o,$o|rv^K^)V'l(^')^|o,^o) 

(0,$o|5|0,$o) 



(2.5) 



(2.6) 



where |0,$o} = |0}|*I'o) represents the noninteracting free Ferm i se a with no p robe particle. Noninteracting electron 
and probe-particle propagators are easily found from Eqs. ( |2.5| ) and (2.6) to be given by the following simple 
expressions: 



and 



G\x,x')^-i ($olTV'(x)V^(x')l$o) 



D°(X,X') = -i (0|TV'(X)V^(^') |0) 



(2.7) 



(2.8) 



respectively. 

We note that the probe-particle propagator is a retarded function, i.e., it is different from zero only if i > t'. As 
a conse quen ce, probe-particle bubbles do no contribute to the diagrammatic expansion. Therefore, the expansion 
of Eq. ( |2.6| ) does not depend on whether the prob e pa rticle is a fermion or a boson, and there is no probe-particle 
contribution to the denominator of Eqs. ( ^ ) and (2.6). 
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A. Scattering approach 



Let us consider the process corresponding to the creation of a single electron- hole pair, where the system is carried 
from an initial state Aj |0, $o) to a final state Oj^ai-^A'j |0, $0)- The scattering-matrix element for this process is 

_ {0,<i>o\af,alAfSAl\0,<i>o) 

(0,$o|^|0,$o) • ^'-^^ 

Similarly, one may consider a double excitation, in which the system is carried from an initial state AI |0, $0) to a 
final state a^j:_^aj^ai^ai2A'j |0, $o)- The matrix element for this process is 

^ _ (0,'i>o|«/, Q/.< alAjSAl\0,'i>o) 



After i ntro ductio n of the Hamiltonian of Eq. (2.1) into Eq. (2.4), the matrix elements 5'/,/i;i,ii and S fj-^j^-AAiM 
of Eqs. ( pj| ) and (2.10) can be expanded in powers of the coupling constant e^. Then, the use of Wick's theorem 
yields explicit expressions for the various contributions to this expansion, in terms of the noninteracting propagators 
G°(a;,x') and D''\X, X'). Introducing standard Fourier representations and taking the free-particle states to be 
momentum eigenfunctions, all contributions to the scattering-matrix can be derived from scattering-like Feynman 
diagrams, as follows: 

1. Draw all distinct scattering diagrams, in momentum space. All particle lines must be directed. Different ways of 
directing them that are not topologically equivalent give distinct contributions. Exclude probe-particle bubbles. 

2. Assign momentum and energy to all particle and interaction lines, so that the sum of the four-momenta entering 
a vertex equals the sum of four-momenta leaving the vertex. 

3. Include an overall factor 27r (5k (5(fc*'), which represents total momentum and energy conservation. 6]^ is the 
Kronecker 6 symbol and S{k'^) is the Dirac S function. 

4. For every external particle line include a factor V^^^^. 

5. For every internal electron line include i G^, where is the noninteracting one-electron propagator in momentum 
space: 

(k, fc*^) being the four-momentum of the particle, rik the occupation number, [nk = Q{qF — IqI), where qp is the Fermi 
momentum], and cjk = k^/2. 

6. For every internal probe-particle line include a factor i D^, where is the noninteracting probe-particle propagator 
in momentum space: 

D° = (2.12) 

(p,p°) being the four-momentum of the particle, and usp — p'^/(2Af). 

7. For every probe_particle-electron and electron-electron interaction line include a factor i Z\ Wq and —i Wq respec- 
tively, Wq being the Fourier transform of the bare Coulomb potential. 

8. For every electron loop include a factor —2. 

9. Integrate over free four-momenta, / ci^(7/(27r)'*. 

Since all scattering-matrix elements include delta functions accounting for momentum and energy conservation, one 
may factorize them as follows 

Sfji-i,ii = 27r5p^_p,_k/^+kii f^i^Pf- ^Pi - ^kf^ + uJk.jTfj^.i^i^ (2.13) 

and 

Sfji,f2;i,ii,i2 = 27r 5p^._p,_k^'^-k,2+k,i+ki2 S{u!p^ - UJp^ - UJ]^^^ - [Jk^^ + ^^'l + ^^'2 ) ^/ Ji J2 ,i2 i (2-14) 

where i^ii,i2jij2 ^^.d pij represent the initial and final momenta of target electrons and probe particle, respectively, 
with energies t^k......^..^, = and Wp,^ = p2^/(2M). 
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The probabilities '^^'"sio g^j^^j ^double ^^j. ^j^g probe particle to transfer four-momentum q {q^ > 0) to a FEG by 
creating single and double excitations are derived by summing the matrix elements over all available initial and final 
electron states and all final probe-particle stateaa. 

^single ^ 4^ ^ „^ (1 _ nu+q) |rq,k(ft)l' S (g" + c^k - c^k+q) <5 [g" - q • V + <zV(2M)] (2.15) 

k 

and 

^double = StT ^ /dgj ^ ^ nk, (1 - nk,+qjnk, (1 - ^k.+q-qj |Tq,q„k„k.(p»)P 
qi ki k2 

(gj + ^k, - ^k,+qO S (qO - + ^k. - ^k.+q-qj (5 [q" - q • V + q^ i2M)] , (2.16) 

where v represents the velocity of the probe particle. 

In these equations recoil has not been neglected. Moreover the quantum character of the probe particle is implicit 
in the T-matrix elements, which include the probe-particle propagator D^(p). Therefore they generalize the results 
of Ref. H to the case of an arbitrary distinguishable probe particle. 

Hence, the total decay rate of the probe charge is given by the following expression: 

/>OC 

r-\p) = E / hr^'^ + 7^'"= + •••] • (2.17) 

q 

The average energy lost per unit length traveled by the pro be pa rticle, i.e., the so-called stopping power of the 
target is obtained by inserting q^ /v inside the integrand in Eq. ( ^.17 ): 



It is well known that the decay rate and the energy loss cannot be computed by simply evaluating the lowest- 
order tree-level Feynman diagrams, because of severe infrared divergences due to the long-range Coulomb interaction. 
Instead, one needs to resum electron-loop corrections and expand the scattering matrix in terms of the dynamically 
screened Coulomb interaction. 

Direct, linear, and quadratic contributions to the screened interaction are represented in Fig. 1. Dashed lines 
[— iwq] represent the bare Coulomb interaction. The full bubble and triangle, denoted ixq and —2Yq-^^q^^ represent 
the sum of all possible Feynman diagrams joining two and three points, and thus correspond to the Fourier transform 
of time-ordered density correlation functions of the interacting FEG : 

Xg= y"d4xie"^['i-('^^-"-^)-«"(*^-*^)]x(xi,X2) (2.19) 

and 

y,,,,, = j d^xi J d^x2e-'i'^'-^'''-'''^-^"^''~'''^]e~'i^'''+'''^-^^^^^ (2.20) 

with 



X{x,x') = -i (^'o|TpH(a;)pH(a;')|*o) 



(2.21) 



and 



Y{x,x',x") = --{■fo\TpH{x) ph{x')ph{x")\^o)- 



(2.22) 

Ph (x) — n is the exact Heisenberg 



Here, |\l/o) represents the normalized exact many-electron ground state, and ph{x 

electron-density fluctuation operator, both in the absence of probe particle. 

Introducing complete sets of energy eigenstates between the Heisenberg fields of Eqs. (2.21) and ( 2.22| ), one obtains 
spectral representations for Xq a-nd Yq-^^^q^. We find: 
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1 



1 

- ^nO + i Vq° 



(2.23) 



and 



Y — 



(Pqi)0n(Pq3)ni(Pq2)m 



(Pq2)0n(Pqi)rii(Pq3)z0 



(q? - w„o + if/go ) {q^ +uJio + i r]gO ) 

(^3 - ^nQ + i?7gO )(q? + + i ?7q? ) 



(4 
(92 



^ 93) 



(2.24) 



where ry, = 77Sgn(g°), = — i?;, (73 = —((71 + (72) and (pq)„; is the matrix element of the Fourier transform of 
the electron-density operator, taken between exact many-electron states of energy En and Ei. 

In the RPA, density correlation functions are obtained by summing over all ring-like diagrams, as shown in Fig. 2, 
thereby neglecting all self-energy, vertex, and vertex-ladder insertions [see Fig. 3]. Hence, 



RPA 



4 



and 



B.PA 



91:92 



91 91.92 —92 —93 ' 



(2.25) 



(2.26) 



where and Y^_^ are the non-interacting FEG density correlation functions and Kq is the so-called inverse dielectric 
function: 



K, 



RPA 



1 



Xq Wq- 

Improvements on the RPA are typically carried out by introducing an effective e-e interactior 

Va 



(2.27) 



.g-Vq(l-Gg), (2.28) 

where Gq is the so-called local-field factor, first introduced by llubbard,0 accounting for all self-energy, vertex, and 
vertex-ladder insertions not present in the RPA. Accordingly, the density correlation ftuactions Xq and Yq-^^^q^ are found 
to be of the RPA form, but with all e-e bare-Coulomb interactions Wq replaced by Wglla: 



and 



Xq = xl+X^VqXq, 



K Y^ k k 

91 -'91,92 -'^-92 -f^-93' 



where the so-called test_charge-electron inverse dielectric function has been introduced 



00 



1 



Xq Vq. 



(2.29) 
(2.30) 
(2.31) 



This inverse dielectric function screens the potential generated by a distinguishable te st cha rge and 'felt' by an electron, 
whereas the so-called test_charge-test_charge inverse dielectric function Kq of Eq. (2.27) screens the potential both 
generated and 'felt' by a distinguishable test charge: 



Kq = 1 



Xq Vq- 



(2.32) 



Now we proceed to expand the matrix elements of Eqs. ( |2.9| ) and ( 2.10 ) in powers of the dynamically screened 
Coulomb interaction. At this point, we will only introduce self-energy and vertex insertions that can be described 
with the use of a static local-field factor [Gq — > Gq,o]- Within this approximation, all processes corresponding to 
the creation of single and double excitations can be represented by diagrams of_Fig. 4. [as in the RPA there are no 
contributions, up to second order in Zi, from triple and higher-order excitationsO] . Thus, one finds: 



T^M = iZiV-^v^Kq + ZlV-^ 



(2^) 



4 [2 Vq Wqj Wq-qi ^p-qi ^q.-qi 



+Vqi Kq^ "q-qi ^q-qi ^p-qi iP^- 



-91 



G°k+q-qi) 



(2.33) 
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and 





k2+q 



T'q.qi.ki.ks ^ iZlV ^ [2 Wq Vq^ Vq-q^ Yq-qi 

+Vq Kq Wq-qi -^^q-gi (Cfci+g + Gl^_q^qJ + Vq^ Kq Uq^ Kq^ {G\ 
— iZlV ^ Uqj Kq^ ^q-qi -^<?-gi -^'p-g^ • 

As static local-field factors arc known to be real [ImGq_o = 0], one easily finds 

IvciKq = Wq |^g|^Inix° 

and 

lvakq^iq\kq\HmX% 



(2.34) 



(2.35) 



(2.36) 



where Kq and Kq are the inverse diel ectri c functions of Eqs. (2.27) and ( 2.32| ), with the density correlation function 

Xq being given in both cases by Eq. ( 2.29 ). 

Introduction of Eqs. ( |2.33| ) and ( p. 34 ) into Eqs. (2.15) and (2.16) yields the following results, valid up to third 
order in the probe_particle-electron screened interaction; 



single 
'9 



2ZlV'^Vq{lmKq + AZi 



(2^) 



4 ^qi ^q-qi 



lmK„ 



Im ( K* Kq, Kq_q, D 



P-qi ^q-.qi j\ j 



X(5[(7"-q-v + gV(2Af)]e(<z°) 



and 



double 



16 Zf vq J ^ «q, i;q_q, { luii^,, ImKq^q, Re (X, D";_^, Yl 



-linKq^q, Re 



K K* 



{Df_q, + + Iq,,q] } 5 - Q • V + / {2M)\ 



xe(gj)e(g"-<z?), 



(2.37) 



(2.38) 



where we have defined the function /„ „, as 



\-^q.,qi ^" -^q,q—Qi ^ {^q,qi ~^ '^f 



q-.q-qi Ji ' 



(2.39) 



with 



-2^y-i^ nk(l-nk+q) 



S (g° + t^k - ^k+q) S (q° - CJk + t^k+q) 



9l + ^k — ^^k+qi 



(2.40) 



and 



Jq^q, = 2tT^ V ^ nk (1 - nk+q)(5 + CJk - t^k+q) 

k 

X [(1 - '^k+qi) <5 {Qi + t^k - ^k+qi) - '^k+q-qi S {q° - Qi + UJk - ^^k+q-qi)] 



(2.41) 



Introduction of Eqs. (2.37) and ( ^.38| ) into Eqs. (2.17) and ( |2.18| ) yields the total decay rate and the average 
energy loss of arbitrary particles that are distinguishable from the electrons in the Fermi gas, with inclusion of static 
many-body local-field effects. 



In order to compare our result with previous work, we consider now the case where the probe particle is very heavy 



(M >> 1) and recoil can be neglected, i.e., Wp 



q • V. It this approximation, the principal part of the 



noninteracting probe-particle propagator is found to give no contribution to the integrals of Eqs. (2.37) and (2.38), 
and one finds 



6 



ImK, 



K K 



Re 



K K T 

J-^qi ^^q-qi ^q,qi 



} <5(g"-q.v)e(0 



(2.42) 



and 



double 

iq 



{2tt 



4 ''qi i^q-qi 



91 



-2 ImK 



<5 (g? - qi • v) <5 (gO - q . v) e(g?) e(g° - q°). 



(2.43) 



These decay probabilities, which account for the existence of many-body static local-field corrections, coincide in 
the RPA [Gq = 0] with those derived in Ref. ^ by treating the probe particle as a prescribed source of energy and 
momentum. 

Simplified expressions for the total decay rate and the average energy loss of heavy (M >> f ) probe particles can 
be obtained with the aid of the following relationship, obtained in Ref. ^, which relates the imaginary part of the 
noninteracting density correlation function Y^^_^^ with the function Hq^q^ of Eq. ( p. 40 ): 



^-92,91 + ^-93:92 + — > 93)] 



(2.44) 



Introduction of Eqs. (2.42) and (2.43) into Eqs. (2.17) and (2.18) yields, after some algebra, 

i;q I -ImKq + 47r Zl J j—^ Vq, Wq_qi S - qi • v) 
X [fl{q, qi) + f2{q, qi) + f^iq, qi) + fl{q, qi) + fiiq, gi)] } <5 (g° - q • v) 6(9") 



(2.45) 



an< 



dE __4Tr 



where 



and 



Itt 2 /" d q of f d*qi r / \ 

V^^J (2^^ Vqi-lmKq + 4nZi J Vc^, Vc^-c^, S [q^ - qi ■ v) 

X [/l(g,9l) + ./2(g, qi) + n{q,qi) + h{q, 9l)]} ^ (g° - q ' v) e(g"). 



/!(<?, 9i) = Irakq ReF,°_,^ YleKq^ ReKq_q^, 

f2{q, qi) = Rekq Hq,q^ ReKq^ RoKq^q^ , 
fsiq, 91) = -2 Imif, Hq.^q ImKq, ReKq.q, , 
/a (9, = -Rei^g iJq.qi ImKq^ ImKq-q^ , 

/4('?, gi) = - o Im-^9 Im^qi Im^^^^, 



Mq,qi) = Im{KqK*K, 



'■q ^^qi ^^q-qi I ''9-91,-91 • 



(2.46) 

(2.47) 
(2.48) 
(2.49) 
(2.50) 

(2.51) 
(2.52) 



Withi n RP A, the inverse dielectric functions Kq and Kq coincide, and Eq. (2.46) reduces to the result of Ref. |[ 

Eq. 2.45 has not been reported before, even within the RPA approximation. In next section we will show that it is 
equivalent to the result reported in Ref. |l^, where a derivation of the decay rate as the imaginary part of the on-shell 
self-energy of the probe particle was sketched briefly. 



B. Self-energy approach 

Since we are considering the interaction of a moving probe particle with a spatially uniform electron gas, invariapt 
under translations, the exact probe-particle propagator can be written in the form of an algebraic Dyson's equationtil 
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Dn 



V 



Dp Tip Dp , 



(2.53) 

which defines the self-energy Ep of the probe particle. With the aid of Eq. ( ^.12| ), Dyson's equation can be solved 
explicitly as 



Dp = 



(2.54) 



The energy and lifetime of the excited state (quasiparticle) obtained by adding a particle to an interacting ground 
state are determined by the poles of the analytical continuation of the one-particle Green function. Hence, the energy 
of the quasiparticle is Wp -I- ReEp, and the probability for it to occupy a given excited state decays exponentially in 
time with the decay constant 



-2ImE, 



(2.55) 



with the self-energy calculated at the pole of the one-particle propagator Dp. 

The self-energy can be represented diagrammatically as the sum of the so-called proper self-energy insertions, i.e., 
all Feynman diagrams that cannot be separated into two pieces by cutting a single particle line. Since the probe 
particle, of charge Zi, is assumed to be distinguishable from the electrons in the Fermi sea, the self-energy may be 
expanded in powers of Zi, diagrams of order containing n — 1 probe-particle propagators. For a homogeneous 
electron gas, contributions from the uniform positive background are canceled by the sum of the so-called 'tadpole' 
diagrams; therefore, after resuming all electron-loop corrections, the self-energy of the probe particle can be represented 
diagrammatically up to third order in Zi as in Fig. 5. The sum of the first two diagrams represents the so-called GW 
approximation, and the third diagram accounts for Zf corrections to the decay rate of the quasiparticle. One finds: 



Ep — z Z-^ 



(27r)4 '^'1^^'-' 



(1 + Xq Vq) -2iZl 



Al 

(2^; 



4 Dp-qi Dp-qj^q^ Yq^-q-^ 



'qi i^q-qi 



(2.56) 



where Xq a nd Y"q , represent the exact density correlation functions of the interacting FEG, as obtained from Eqs. 
( ^.231) and ( |2.24)| ), respectively. 



If the probe particle is an ion (M >> 1), the propagator Dp and the energy p° entering Eq. (2.56) can be safely 
approximated by the noninteracting propagator Dp and energy uj-p. Furthermore, recoil can be neglected, and one 
easily finds 



D^ 

p-q 



q • V — 2 ?; 



(2.57) 



In order to exploit the symmetry properties of Yq-^ ^q^ it is useful to rewrite the retarded probe-particle propagator in 
terms of its Feynman version as follows, 



D°=- 



1 



q • V 



''Vq 



-2z7r(5(g°-q-v)e(q°). 



(2.58) 



Introducing Eq. (2.58) into Eq. ( ^.56 ) and noting that the time-ordered density correlation functions Xq ^-nd ^gi,g2 
are invariant under the changes {q'^ — > —q^) and {qi — > —qi, ^2 ~* ~92)j respectively, some work of rearrangement 
leads us to the following expression: 



= 2 TT 



dq^ 



4 



(1 



Xq Vq)-^TrZl J Yq^^q, Vq^l^q-qi 6{q^ - qi • v) 



X(5(g"-q-v)e(g°). 



(2.59) 



Within RPA, the density correla tion functio ns X a and Yq-^^q^ are those given by Eqs. (2.25) and (2.26). Beyond 



RPA, they are obtained from Eqs. ( 2.2£ ) and ( 2.3Cj), in terms of the noninteracting density correlation fimctions [Xq 
and ^J^^gj] and the effective e-e interaction of Eq. ( 2.28| ). Hence, introduction of Eq. ( 2.59| ) into Eq. ( 2.55 ) yields the 
following expression for the decay rate: 



n Z\ 



(2^) 



^\^[KqYl_^^kq,kq^q, 



fqiWq-qi <5(9? " Ql ' v) 



(2.60) 
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where Kq and Kq represent the inverse dielec tric functions of Eqs. (2.27) and ( ^.32 ), with the density correlation 
function Xq being given in both cases by Eq. ( 2.29 ). 

The equivalence of Eqs. ( 2.45| ) and (2.6C) follows from the expansion of the imaginary part oi KqYq _q_^ Kq^Kq^q^ 

m Eq. ( ^!60| ), and the use of Eq. ( 2.44[ ) and the symmetry properties of Yq^^q^ and K{q^,(\). However, while ( |2.45| ) 
has been derived by only introducing self-energy and vertex inse rtion s that can be d escribed with the use of a static 
local-field factor, we have now demonstrated that either Eq. ( 2.45| ) or Eq. ( 2.60| ) can be used with inclusion of 
many-body dynamic local-field corrections. 

Finally, we note that although both Eq. ( |2.45| ) [derived from Eq. ( |2.17D ] and Eq. ( |2.60D [derived from Eq. ( ^.551 )] 
represent the total decay rate, the integrands of these integral representations do not necessarily coincide with the 
probability [^^™ff'^ -|- ^double _|_ _j fQj- ^j^g probe particle to transfer four- momentum q to the FEG. Consequently, the 
stopping power of th e FE G for the probe particle [see Eq. ( ^.46 )] cannot be obtained by simply inserting q'^ /v inside 
the integral of Eq. ( 2.45| ) or Eq. (2.6C), and the knowledge of the self-energy alone is not, therefore, sufficient to 
calculate the stopping power. 



C. Quadratic response 



In Ref. 10 the stopping power of a heavy probe-particle was calculated using the framework of quadratic response 
theory. It was found that 



dE 
dx 



47rZ?/^<Z%,5(g"-q-v) e{q") 

Imi^f + 27r Zi 1 ^ Im (xf X« v^, v^^^, S (g? - qi ■ v) 



(2.61) 



where K^, Kq and Y^^q^ represent the retarded counterparts of Kq, Kq and i^gi, ga- 
in order to compare this result to the one quoted in Sec. 2 we must recall the relationship between the time-ordered 
and retarded functions. 

In our case of interest, differences in the inverse dielectric functions arises from differences in the FEG density 
correlation fimction x". We have: 



Rexf° = Rex°, 



(2.62) 



Imxf ° = sgn((j°)Imx°, 



The relationship between Yq-^^q^ and Yj^ is best analyzec 



(2.63) 



using their spectral representations. Y^^ q^ has the same 
structure as Eq. ( 2.24 ), with rjq replaced by a positive rj^ This property leads to the following relations between 
imaginary and real part of the time-ordered and retarded Y^^ ^2 functions: 



Re 



i?,,0 



J-q2.q3 + J-q3,q2 j 



(2.64) 



and 



= 2 [sgn('Z?)^«.,-g. -sgn(<Z2°)i?-, 



Sgn((7§)i?-g3,g2 



(92 



93 J 



(2.65) 



After some work of rearrangement, and taking into account the symmetry properties of the functions involved, we 
find that E q. (2.61 ) coincides exactly with Eq. ( 2.46| ). As in the case of the decay rate of Eq. ( 2.6C ), we find that 
both Eqs. (2.46) and (2.61) can be use d wit h inclusion of many-body dynamic local-field corrections. We also note 
that within RPA both Eqs. ( |2.46| ) and ( |2.6l| ) reduce to the result derived in Refs. 



III. CONCLUSIONS 



We have developed various many-body theoretical approaches to the quadratic decay rate and energy loss of charged 
particles moving in an electrons gas, with inclusion of short-range XC effects. 
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We have carried out a perturbative formulation of the scattering matrix to derive general expressions for both the 
total decay rate and the average energy loss of arbitrary moving charged particles that are distinguishable from the 
electrons in the Fermi gas. Simplified expressions for these quantities have been obtained in the case of heavy probe 
particles [M » 1). The total decay rate of heavy particles has then been rederived from the knowledge of the 
probe-particle self-energy and it has been proved that the stopping power of the heavy particle agrees with the result 
deduced using quadratic response theory. Comparison of the different formalisms for a heavy particle suggests that 
our results in the scattering formalism can be used with full inclusion of many-body dynamic local- field corrections. 

It has also been shown that while the first-order contributions to the energy loss may be obtained from the total 
decay rate by simply inserting the energy transfer inside the integrand of this quantity, this procedure cannot be 
generalized to the description of the second-order energy loss. Since response theory is only valid for heavy particles, 
this implies that the stopping power of light particles must be calculated using scattering theory. 
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FIG. 1. (a) Direct, (b) linear, and (c) quadratic contributions to the screened interaction. Dashed lines represent the bare 
Coulomb interaction, —i «q. Two- and three-point loops represent time-ordered density correlation functions, i Xq and —2 Yqj , 
respectively. 
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FIG. 2. (a) The interacting RPA two-point density correlation function, represented by a full bubble, is obtained by summing 
over the infinite set of diagrams that contain a string of empty bubbles, (b) The interacting RPA three-point density correlation 
function, represented by a full triangle, is obtained by summing over the infinite set of diagrams that combine three strings of 
empty bubbles through an empty triangle. 



FIG. 3. (a) Self-energy, (b) vertex, and (c) ladder insertions, which are neglected within RPA. 



FIG. 4. Diagrammatic representation of the scattering-matrix elements of (a) Eq. (2.33) and (b) Eq. (2.34). Thick and 
thin solid lines represent noninteracting probe-particle and electron propagators, i_Dp and iG^l, respectively. Wiggly lines 
represent the screened e-e Coulomb interaction, ~i Vq Kq , and thick discontinuous lines represent the screened interaction 
between electrons and probe particle, —iVqKq 



FIG. 5. The probe-particle self-energy, up to third order in Zi. Thick solid lines represent the exact probe-particle propagator, 
iDp. Dashed lines represent the bare Coulomb interaction, — z«q. Two- and three-point loops represent time-ordered density 
correlation functions, iXg and — 2Fqj,52, respectively. 
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Figure 5 




